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We set up a Chern-Simons theory for the entanglement of 
two polymers Pi and Pi , and calculate the second topological 
moment (m 2 ), where m is the linking number. The result 
approximately to a polymer in an ensemble of many others, 
which are considered as a single very long effective polymer. 



I. THE PROBLEM 

Consider two polymers P\ and P2 which statistically 
can be linked with each other any number of times m — 
0, 1, 2, . . . . The situation is illustrated in Fig. |] for m = 
2. 



(m 2 ) = 




FIG. 1. Closed polymers Pi, P2 with trajectories C\, C2 re- 
spectively. 

We would like to find the probability distribution of 
the linking numbers m as a function of the lengths of Pi 
and Pj. As a first contribution to solving this problem 
we calculate, in this note, the second moment of this 
distribution, (to 2 ). The self-entanglements are ignored. 

The solution of this two-polymer problem may be con- 
sidered as an approximation to a more interesting phys- 
ical problem in which a particular polymer is linked to 
any number N of polymers, which are effectively replaced 
by a single long "effective" polymer jlj] . 

Let G m (xi, X2; Li, L 2 ) be the configurational probabil- 
ity to find the polymer Pi of length L\ with fixed coin- 
ciding end points at xi and the polymer P2 of length L2 
with fixed coinciding end points at X2, entangled with a 
Gaussian linking number to. 



The second moment 
integrals 
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is defined by the ratio of 
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/ d 3 Xid 3 x 2 J_™ dm m 2 G m (xi, x 2 ; L 1} L 2 ) 
J d 3 xid 3 x 2 f_™ dmG m (x 1 ,x 2 ;L 1 ,L 2 ) 



(1) 



performed for either of the two probabilities. The in- 
tegrations in d 3 xid 3 x 2 covers all positions of the end 
points. The denominator plays the role of a partition 
function of the system: 



//•-too 
d 3 x 1 d 3 x 2 / dmG m {x 1 ,x 2 ;L 1 ,L 2 ) (2) 
J — 00 

Due to the translational invariance of the system, the 
probabilities depend only on the differences between the 
end point coordinates: 

G m (xi,x 2 ; Li,L 2 ) = G m (xi - x 2 ;Li,L 2 ) (3) 

Thus, after the shift of variables, the spatial double inte- 
grals in ([!]) can be rewritten as 

J d 3 x 1 d 3 x 2 G m (x 1 -x 2 ;Li,L 2 ) 

= V J d 3 xG m (x;L u L 2 ), (4) 

where V denotes the total volume of the system. 



II. POLYMER FIELD THEORY FOR 
PROBABILITIES 

The linking number for the two polymers is given by 
the Gauss integral 



7 G (Pi,P 2 ) = -L i I [dxi x dx 2 ] X ' " X2 



4tt 



Pi JP: 



xi - x 2 



(5) 



It takes the values to = 0, ±1, ±2, .... With the help of 
two vector potentials Ai and A 2 , the phase factor e lmA 
can be obtained as a result of a local gauge theory of the 
Chern-Simons type: 



-Acs — kJ p dxi-Ai — ^J p dx 2 -A 2 



x e 



(6) 



where -4ns is the action 



1 



Acs = in 



xt 



txvp 



A 1 d u A2, 



(7) 



Indeed, the correlation functions D^{x, x') of the gauge 
fields are 



(A?(x)Al(x>))=0, (^(x)^(x')>=0, 
K(x)^(x')) = 



d 3 P jp(x-x') i e ^yk x 



(2tt) 3 
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(x - x') 

\Z - Y '|3 



(8) 
(9) 

(10) 



such that the functional integral on the right-hand side 
of (||) produces directly the phase factor e l,Ia ^ Pl,p ^ x with 
the eigenvalue e mlX . 

We can select configurations with a certain linking 
number m from all configurations by forming the inte- 
gral d\e~ vmX over this quantity. 

The most efficient way of describing the statistical 
fluctuations of the polymers P\ and Pi is by two com- 
plex polymer fields ^^(xi) and ?/>2 2 ( x 2) with n\ and n% 
replica {a\ = 1, . . . , m, 02 = 1, ... , 712)- At the end we 
shall take 711,7x2 — > to ensure that these fields describe 
only one polymer each For these fields we define an 
auxiliary probability Ga(xi, X2; z) to find the polymer 
Pi with open ends at Xi , x^ and the polymer P2 with 
open ends at x 2 ,x 2 . The double vectors x x = (xi,x^) 
and X2 = (x2,x 2 ) collect initial and final endpoints of 
the two polymers Pi and Pi. The auxiliary probability 
G\ (xi , X2; z) is given by a functional integral 



G\ (xi , X2 ; z) = lim / P(fields) 

ni,ti 2 - >0 J 

x ^ ( Xl )^ 01 (x'i)^ 2 (x 2 )^r 2 (^k- 4 



(ii) 



where 2?(fields) indicates the measure of functional in- 
tegration, and A the action governing the fluctuations. 
It consists of kinetic terms for the fields, a quartic in- 
teraction of the fields to account for the fact that two 
monomers of the polymers cannot occupy the same point, 
the so-called excluded-volume effect, and a Chern-Simons 
field to describe the linking number m. Neglecting at first 
the excluded-volume effect and focusing attention on the 
linking problem only, the action reads 



A — Acs + Apoi, 
with a polymer field action 



(12) 



2 

A po i = ]T / d 3 x [|D%| 2 + m 2 |^| 2 ] . (13) 



in which we have omitted a gauge fixing term, which 
enforces the Lorentz gauge. They are coupled to the 
polymer fields by the covariant derivatives 



D i = V + *7 i A i , (14) 

with the coupling constants 71.2 given by 

71 = K 72 = A. (15) 

The square masses of the polymer fields contain the 
chemical potentials zi_i of the polymers: 



mi = 2Mzi. 



(16) 



They are conjugate variables to the length parameters L\ 
and I/2, respectively. The symbols vf^ collect the replica 
of the two polymer fields 

= *< = (#V ( 17 ) 

and their absolute squares contain the sums over the 
replica 



|D%| 2 = J2 |Dty"f . I^l 2 = E l^l 2 - (18) 

Having specified the fields, we can now write down the 
measure of functional integration in Eq. ( |Tl| ) : 

©(fields) = J VA^VA^V^iV^*iV^ 2 V%. (19) 

By Eq. (JsJ) , the parameter A is conjugate to the link- 
ing number m. We can therefore calculate the de- 
sired probability G m (xi, X2; L%, L2) from the auxiliary 
one Ga (xi , X2 ; ,?) by the following Laplace integrals z — 
(zi,zi): 

G m (x 1 ,x 2 ;L 1 ,L 2 )= Hm / d Jl^l e i^ 1+Z2 L 2 ) 
"i-"! Jc-ioo 2m 2m 

/OO 
d\e~ imX G\ (xi, x 2 ; z) . (20) 
-OO 



III. CALCULATING THE PARTITION 
FUNCTION 

Let us use the polymer field theory to calculate the 
partition function (||). By Eq. (pp|), it is given by the 
integral over the auxiliary probabilities 



Z = / cPxicPx2 lim 



c+00 



2m 2m 



/+00 
dm d\e- lmX Gx(xi,x 2 ;z). (21) 
-00 J —oo 

The integration over dm is trivial and gives 2tt5(X), en- 
forcing A = 0, so that 



2 



z 



d xicrx2 Urn 



c+zoo 



dz\dz% 
2iri 



(z 1 L 1 +z 2 L 2 ) 



x G A =o (xi,x 2 ;z) 



(22) 



To compute Ga=o (x"i , x 2 ; z) we observe that the action 
A in Eq. (O) is quadratic in A. Let us expand A as 



.4 = A + XAi + X A 2 



(23) 



where 
A = Acs 



dh 



|Di^i| 2 



|V* 2 | 2 
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a linear coefficient 

At = J d 3 x j 2 (x) • A 2 (x) 

with a "current" of the second polymer field 

j 2 (x)=**5(x)V* 2 (x), 
and a quadratic coefficient 

A,= \J d 3 x A 2 |* 2 (x)| 2 . 

With these definitions we can rewrite (|23| ) as 

G A=0 (x 1; x 2 ;z)= /" P(fields)e-- 4 ° 

x^ 1 ai (xi)V'r i (x 1 )V' 2 2 (x 2 )v 2 2 (x') 



(24) 



(25) 



(26) 



(27) 



(28) 



From Eq. (£4|) it is clear that G\=o (xi , x 2 ; z) is the prod- 
uct of the configurational probabilites of two free poly- 
mers. 

Note that the fields \E , 2 ,\E , 2 are free, whereas the 
fields ^i, are apparently not free because of the cou- 
plings with the Chern-Simons fields through the covari- 
ant derivative D 1 . This is, however, an illusion: the fields 
A 1 ^ have a vanishing diagonal propagators (A* A\) = 0. 
integrating out A% in (p8|), we find the flatness condition: 

e^ p d v A\ t = 0. (29) 

On a flat space with vanishing boundary conditions at 
infinity this implies A^ = 0. As a consequence, the func- 
tional integral ((H) factorizes 

G A =o (xl,x 2 ; z) = G (xi - x' x ; zi)G (x 2 - x' 2 ; z 2 ), (30) 

where Go(x^ — xj; z,) are the free correlation functions of 
the polymer fields: 

Go(x t - x' i5 zi) = (^(x^Vr^x-)}- (31) 

In momentum space, the correlation functions are 



(^(kSr(K)) = 5 (3) (ki + kj) (32) 

K -t" TO„- 



such that 



G (x 4 -x, t ;z. i ) = J — 



d 3 k e ik- x 1 



3 k 2 +TO 2 ' 



(33) 



and 



G (x 4 - x' 4 ; ^) = / -^e 2 - i -G (x J - x^) 



(34) 



4V2M \2irJ 
Thus we obtain for (p2]): 

Z = 2n J d 3 x 1 d 3 x 2 

x lim G (xi - x' 1 ;ii)G (x 2 - x' 2 ;L 2 ) (35) 



The integrals at coinciding end points can easily be per- 
formed, and we find 



(36) 



It is important to realize that in Eq. ([H|), the lim- 
its of coinciding end points x£ — > Xf and the inverse 
Laplace transformations do not commute unless a proper 
renormalization scheme is chosen to eliminate the diver- 
gences caused by the insertion of the composite operators 
IV'MI 2 - This can be seen for a single polymer P. If we 
were to commuting the limit of coinciding end points with 
the Laplace transform, we would obtain 



/ e zL lim G (x-x';z) 

Jc-ioo 27T x'^x 



c+zoo 



dz 



2m 



T e zL G (0,z), 



(37) 



where 

G (0;z) = (|^(x)| 2 ). (38) 
This expectation value, however, is linearly divergent: 

d 3 k 



(l^(Xa)| 2 ) = 



(39) 



IV. CALCULATION OF NUMERATOR IN 
SECOND MOMENT 

Let us now turn to the numerator in Eq. (|l|): 
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N = /"d 2 xid 3 x 2 / dm m 2 G m (xi, x 2 ; Li, L 2 ) • (40) 

We set up a functional integral for N in terms of the aux- 
iliary probability Ga=o (x"i, x 2 ; z) analogous to Eq. (fj"l|): 



TV = / d 3 xid 3 r 2 



<im m 2 lim 



d^j dz2 
2m' 2ni 



= (ziil+Z2-L2) 



-imX 



(41) 



The integration in dm is easily performed after noting 
that 



dm m 2 e~ unX G x (xi,x 2 ;z) 

i 



(42) 



After a double integration by parts in A, and an integra- 
tion in m. we obtain 



N 



/d 3 x x d 3 x 2 lim (-1) f M °^£i^i e ( Zl ii+^ 2 
7 i 7 c -ioo 2wi2m 



dA 5(A) 



dA 2 



G A (xi,x 2 ; f) 



(43) 



Performing the now trivial integration in dX yields 



N - 



c-Hoo 



'/'x.r/'x, lim (-1) / ^i^e 
' 2m 2m 



{Z1L1+Z2L2) 



d 2 

-gj^Gx (xi,x 2 ; z) 



(44) 



A=0 



To compute the term in brackets, we use again Eq. ( pa ) 
and Eqs. @ and find 



N = / cf xi<rx 2 lim 



c+ioo 



^£l^f2 ( Z1 L 1+Z 2L 2 ) 



27TI 27TZ 



x / P(fields) exp(-A)l^ 1 (xi)| 2 |V 2 2 (x 2 )|^ 

yd 3 xA 2 -**v* 2 ^ +iy"d 3 xA 2 i* 2 



and interacting parts 

^° = y" d 3 xj 1 (x)-A 1 (x) 

with a "current" of the first polymer field 
jx(x) = i*J(x)V*i(x), 

and 



(47) 



(48) 



(49) 



Expanding the exponential 

A _ X!+-4o+-4o = P -4o 



1 , K) 5 



(50) 



and keeping only the relevant terms, the functional inte- 
gral ( |45| ) can be rewritten as a purely Gaussian expecta- 
tion value 



N = K< 



d 3 xi<i 3 x 2 lim 



c+ioo j j 

C(Zl Q^2 c ( zlLl +z 2 L2) 

_,•„„ 27TI 27TI 



x / ©(fields) exp(-^)|^ 1 (xi)nV' 2 2 (x 2 )| 



|2|*A°2 f 



)2 
+iy > d 3 xA 2 i* 1 | 



d 3 xA 2 • **V* 2 ) + J / d 3 x A 2 |* 2 



(51) 



Only four diagrams shown in Fig. 
Eq. (fjj). 



contribute in 



FIG. 2. Four diagrams contributing in Eq. (|5l|). The lines 
indicate correlation functions of ^-fields. The crossed circles 
with label i denote the insertion of ^(x;)! 2 . 



(45) 



Note that the initially asymmetric treatment of poly- 
mers Pi und P 2 in the action (jl^) has led to a completely 
In this equation we have taken the limits of coinciding symmetric expression for the second moment, 
endpoint inside the Laplace integral over £i,z 2 . This Only the first diagram in Fig. || is divergent due to the 

will be justified later on the grounds that the potentially divergence of the loop formed by two vector correlation 
dangerous Feynman diagrams containing the insertions functions. This infinity may be absorbed in the four-* in- 



of operations like |*j| vanish in the limit n\, n 2 — > 0. 



teraction accounting for the excluded volume effect which 



In order to calculate (Ha), we decompose the action we do not consider at the moment. No divergence arises 



into a free part 



A = Acs 



d 3 x 



ID 1 *!! 2 + 1 v* 2 | 2 



from the insertiona of the composite fields |*i(xi)| 2 . In 
this respect, the disconnected diagrams shown in Fig. || 
are potentally dangerous. 



(46) 



4 



FIG. 3. Four diagrams contributing in Eq. (|5l|). The lines 
indicate correlation functions of ^-fields. The crossed circles 
with label i denote the insertion of |^i(xi)| 2 . 

But these vanish in the limit of zero replica indices 
ni,m — ► 0. 



V. CALCULATION OF FIRST FEYNMAN 
DIAGRAMS IN FIG. || 

From Eq. @ one has to evaluate the following integral 

AT 1= ^lim f C+t0 ° ^l^l e (^ 1+Z2 L 2) 
4 2m 2m 



J d 3 xid 3 x 2 J (Px'^x'z 



(52) 



K ai (xi)| 2 |^ 



! (x 2 )| 2 (l*i| 2 A?) 



(|^ 2 | 2 Al) 



There is an ultraviolet-divergent contribution which 
should be properly regularized. The system has, of 
course, a microscopic scale, which is the size of the mo- 
momers. This, however, is not the appropriate short- 
distance scale to be uses here. The model treats the 
polymers as random chains. However, the momomers of 
a polymer in the laboratory are usually not freely mov- 
able, so that polymers have a certain stiffness. This gives 
rise to a certain persistence length £o over which a poly- 
mer is stiff. This length scale is increased to £ > £o by the 
excluded- volume effects. This is the length scale which 
should be used as a proper physical short-distance cut- 
off. We may impose this cutoff by imagining the model 
as being defined on a simple cubic lattice of spacing £. 
This would, of course, make analytical calculations quite 
difficult. Still, as we shall see, it is possible to estimate 
the dependence of the integral iVi and the others in the 
physically relevant limit in which the lengths of the poly- 
mers are much larger than the persistence length £. 

An alternative and simpler regularization is based on 
cutting off all ultraviolet-divergent continuum integrals 
at distances smaller than £. 

After such a regularization, the calculation of N\ is 
rather straightforward. Replacing the expectation values 
by the Wick contractions corresponding to the first dia- 
gram in Fig. ||, and performing the integrals as shown in 
the Appendix, we obtain 

V M 2 i 
Ni = — 77r -^{L 1 L 2 )'^ 



4tt (8tt) 6 



x / <fe[(l-s)s]~ 
Jo 



d 3 j 



-Mx 2 /2a(l-s) 



(53) 



1 dt[(l-t)t]-i [ dV- M * 2 /2*(i-*) / d 3 x'[-\ 



The variables x and y have been rescaled with respect 
to the original ones in order to extract the behavior of 
Ni in Li and L 2 . As a consequence, the lattices where x 
and y are defined have now spacings £/^/Li and ti/^/L^ 
respectively. 

The x, y integrals may be explicitly computed in the 
physical limit L\,L 2 » £, in which the above spacings 
become small. Moreover, it is possible to approximate 
the integral in x" with an integral over a continuous vari- 
able p and a cutoff in the ultraviolet region: 



WT J i p 2 



After these approximations, we finally obtain 



8 (8tt) 3 



(54) 



(55) 



VI. CALCULATION OF SECOND AN THIRD 
FEYNMAN DIAGRAMS IN FIG. || 

Here we have to calculate 

N 2 = J lim r +l °° ±l^l e (^L 1+Z2 L 2) 
"IZ° Q Jc-ioo 2?™ 2m 

x J d 3 x x d 3 x 2 J d 3 x' 1 d 3 x'{d 3 x' 2 

x^|^ I (xi)| 2 |V2 2 (x2)| 2 (A 1 .^V^ 1 ) x , 

x(A 1 -^ 1 V^ 1 ) x ,,{A 2 2 \^ 2 \ 2 ) x ,\ (56) 

The above amplitude has no ultraviolet divergence, so 
that no regularization is required. The Wick contractions 
pictured in the second Feynman diagrams of Fig. || lead 
to the integral 



N 2 = -2VL- 1,2 L- l 1 — 



(2tt) 



i ,t 
dt / dt'C{t,t') (57) 
o Jo 



where C(t,t') is a function independent of L\ and L 2 

C(t,t') = [{1 - t)t' {t - t')Y z/2 

x J d 3 xd 3 yd 3 ze- M ^ 2 ' 2 ^ 

5^2, ■ (z + x) - (z + x)^ z 



z P z + x ■ 



(58) 
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As in the previous section, the variables x, y, z have been 
rescaled with respect to the original ones in order to ex- 
tract the behavior in L\. 

If the polymer lengths are much larger than the persis- 
tence length one can ignore the fact that the monomers 
have a finite size and it is possible to compute C(t,t') 
analytically leading to 



^2 = -- 



vl: 



-1/2 



ir 1 



4 3 (2tt) 6 
where K is the constant 

K 



M~ 1/2 V2K, 



(59) 





'3 






'5 


1 




v2' 


\) 














v2' 


2 



X Rf 9 1 



19tt 
384 



0.154, 



(60) 



and B(a,b) = T(a)T(b)/T(a + b) is the Beta function. 
For large L\ — > oo, this diagram gives a negligible con- 
tribution with respect to N±. 

The third diagram in Fig. || give the same as the sec- 
ond, except that L\ and L 2 are interchanged. 



^3 = N 2 \ Ll * 



■L 2 - 



(61) 



VII. CALCULATION OF FOURTH FEYNMAN 
DIAGRAM IN FIG. || 



Here we have the integral 
N 4 = -4k 2 - lim 



2 "i-o 



dzi dz 2 c(ziLi+Z2 l 2) 
2iri 2ni 



x J d^x, f d 3 x[d a x' 2 d 3 x'(d 3 4 (62) 
x^r(x 1 )| 2 |^ 2 (x^)| 2 (A 1 .vI/*Vvl/ 1 ) x; (A 1 .^V* 1 ) x „ 

x (A 2 -*2V* 2 ) X / (A 2 • * 2 V* 2 ) X „ 



which has no ultraviolet divergence. After some effort we 
find 

^4 = ^^ZriL.L,)- 1 / 2 

2-4 6 (2tt) 1iV ' 

x [ ds [ ds' ( dt [ dt'C(s,s',t,t'), (63) 
Jo Jo Jo Jo 

where 

C(s, s'; t, t') = [(1 - s)s'(s - s')}~ 3/2 [(1 - t)t'{t - t')]~ 3/2 
d 3 p 
(2^)3 

f d^x'd^y'e-^^'-^e-^' 2 ' 2 ^ 



pa p/3 pa p/3 

p* P Z P 2 P Z 



x V?„e 



P -Mu' 2 /2t' 



-M(u'-v') 2 /2(t-t') 



(64) 



and x',y' are scaled variables. To take into account the 
finite persitence length, they should be defined on a lat- 
tice with spacing £/y/L~i. Similarly, u', v' should be con- 
sidered on a lattice with spacing ^/y/L^- Without per- 
forming the space integrations d 3 x'd 3 y'<i 3 u'<i 3 v', the be- 
havior of 7V4 as a function of the polymer lengths can be 
easily estimated in the following limits: 
1. L\ > l;Li > L 2 



N 4 oc L^ 1 

2. L 2 > 1;L 2 > Li 

N t oc L- 1 

3. Li,L 2 ^>\, L2/L1 = a = finite 



N 4 (x L x 



3/2 



(65) 
(66) 

(67) 



Moreover, if the lengths of the polymers are consid- 
erably larger than the persistence length, the function 
C(s, s', t, t') can be computed in a closed form: 



Na 



1 



1 



32 



(2tt) 5 (2tt) 3 / 2 V2 



x ds dt(l - s)(l - t)(st) 1/2 
Jo Jo 

x [Li*(l - s) + L 2 (l - t)s}- 1/2 . 



(68) 



It is simple to check that this expression has exactly the 
above behaviors. 



VIII. FINAL RESULT 

Collecting all contributions we obtain the result for the 
second topological moment: 



<m 2 > 



N 1 +N 2 + N 3 + N 4 
Z 



(69) 



with Ni, N A N 3 , 7V 4 , Z given by Eqs. (|36j), (]55j), (|5S|), 
©.and©. 

In all formulas, we have assumed that the volume V 
of the system is much larger than the size of the volume 
occupied by a single polymer, i.e., V S> L\ 

To discuss the physical content of the result (|6^), we 
assume P 2 to be a long effective polymer representing 
all polymers in a uniform solution. We introduce the 
polymer concentration p as the average mass density of 
the polymers per unit volume: 



6 



M 
P= V 

where M is the total mass of the polymers 



M 



N p 

m n — . 

a 



i=l 



(70) 



(71) 



Here m a is the mass of a single monomer of length a, L/. 
is the length of polymer Pk , and N p is the total number 
of polymers. Thus Lk/a is the number of monomers in 
the polymer P k . The polymer P\ is singled out as any 
of the polymers Pk, say F%, of length L\ — Lj.. The 
remaining ones are replaced by a long effective polymer 



P2 of length L 2 
may also write 



a k ^Lk- From the above relations we 



aVp 



(72) 



In this way, the length of the effective molecule P2 is 
expressed in terms of physical parameters, the concen- 
tration of polymers, the monomer length, and the mass 
and volume of the system. Insertin g (|72[) into (|69|) , with 
N u N 2 , 7V 3 , A 4 , Z given by Eqs. IfT©, ©, ©, 
and (|68|). and keeping only the leading terms for V ^> 1, 
we find for the second topological moment (m 2 ) the ap- 
proximation 



2 Ni+N 2 

(m ) « — 
and this has the approximate form 

C x Li KV2L{ /2 



I 2\ a P 

(m) = 



with K of (§3). 



167r i/2 M 2 (27r)4M 3 / 2 



(73) 



(74) 



IX. SUMMARY 

We have set up a topological field theory to describe 
two fluctuating polymers Pi and P 2 , and calculated the 
second topological moment for the linking number m be- 
tween Pi and P 2 . The result is used to calculate the 
second moment for a single polymer with respect to all 
others in a solution of many polymers. 

In forthcoming work we shall calculate the effect of the 
excluded volume. 



X. APPENDIX 

In this appendix we present the computations of the 
amplitudes Ni, . . . , N4. We shall need the following sim- 
ple tensor formulas involving two completely antisym- 
metric tensors e A " /p : 



pa{3 mr/3 _ Sj/3 s:c 



(75) 



(76) 



The Feynman diagrams shown in Fig. g) corresponds to 
a product of four correlation functions Go of Eq. (|3S|), 
which have to be integrated over space and Laplace trans- 
formed. For the latter we make use of the convolution 
property of the integral over two Laplace transforms f(z) 
and g{z) of the functions /, <?: 



c+ioo j pL 

— e* L f(z)g(z) = / dsf(s)g(L - S) (77) 



2ni 



(i 



All spatial integrals are Gaussian of the form 

d 3 xe~ ax+2b *y = (2n) 3 / 2 a - 3 / 2 e b V/a 5 a>Q (78) 



We are now ready to evaluate iVi in Eq. (p2|). Taking the 
limit of vanishing replica indices we find with the help of 
Eqs.©-©: 



p pL\ plj2 p 

N x = J d 3 xi,d 3 x 2 J ds J dtj d 3 x\d 3 x' 2 
x G (xi-x' 1 ; s)G (x' 1 -xx;ii - S) 
x G (x 2 -x 2 ;i)G (x 2 -x 2 ;L 2 -f) — 



Performing the changes of variables 
t' 



s 

Ti 



t 

L~ 2 



XlX x 

x = — — y 



xi — X, s 



x 2 -x 2 



V4 



(79) 



(80) 



and setting x" = x^ — x' 2 , we easily derive (p3|). 

For small £/VTi and ^j^fTa, we d use the approxi- 
mation (|54|), the space integrals can be done the formula 
(f78|). After some calculation one finds the final result of 
Eq. ©. 

The amplitude N 2 calculated quite similarly. Con- 
tracting the fields in Eq. (|5^), and keeping only the con- 
tributions which do not vanish in the limit of zero replica 
indices, we arrive at 

N 2 = J d 3 x x d 3 x 2 J d 3 x' 1 d 3 x l ld 3 x 2 

x / ds ds'G (x' 1 — xi; Li — s) 
Jo Jo 



x V^Go(x! - x'/ ; s')V£, G (x? - x' x ; s - s') 
x D mA (x / 1 -x 2 ) j D 1/A (x' 1 / -x^) 

diG (x 2 - x 2 ; L 2 - i)G (x 2 ~x. 2 ;t) 



(81) 



where L* MI y(x, x') are the correlation functions (p|)-(p^|) 
of the vector potentials. Setting x 2 = ^/T^u + x 2 and 
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supposing that ^/^/L^ is small, the integral over u can 
be easily evaluated with the help of the Gaussian integral 
( [78|). After the substitutions x" = \fL\y + xi x^ = 
\/L\(y — x) + Xi, x' 2 = \/X7(y — x — z) + xi and a 
rescaling of the variables s, s' by a factor , we derive 
Eq. @ with (gjL 

For small £/vLl, the spatial integrals are easily 
evaluated leading to: 

Ar -V2VL- 1/2 L^M-^ 
No = — 

2 (4n) 6 

\ t £ dm - t) J^E*El (82) 

After the change of variable tf — ► t" = t — t 1 , the double 
integral is reduced to a sum of integrals the type 



c(n,m)= / dtt m / dt't' n \ — - — , m. n — integers. 
K ' Jo Jo V 1-f 6 

These can be simplified by replacing i m by dt m+1 /dt(m+ 
1), and doing the integrals by parts. In this way, we end 
up with a linear combination of integrals of the form: 

■„*^T- B (" + 5-5)- (83) 

The calculations of N3 and N4 are very similar, and may 
be omitted here. 
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